Infrared spectroscopy of Landau levels in graphene 
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We report infrared studies of the Landau level (LL) transitions in single layer graphene. Our 
specimens are density tunable and show in situ half-integer quantum Hall plateaus. Infrared trans- 
mission is measured in magnetic fields up to B = 18 T at selected LL fillings. Resonances between 
hole LLs and electron LLs, as well as resonances between hole and electron LLs are resolved. Their 
transition energies are proportional to \/B and the deduced band velocity is c « 1.1 x 10 6 m/s. 
The lack of precise scaling between different LL transitions indicates considerable contributions of 
many-particle effects to the infrared transition energies. 

PACS numbers: 78.67.-n, 76.40.+b, 78.30.-j, 71.70.Di 



Graphene is the newest member in the family of two- 
dimensional (2D) carrier systems, which have shown 
a spectrum of fascinating new physics over the past 
decades. Graphene, a single atomic sheet of graphite, 
represents the ultimate 2D material. Moreover, its elec- 
tronic band structure differs radically from the parabolic 
bands common to all previous 2D systems. In graphene, 
the conduction and valence bands meet at two inequiv- 
alent, charge-neutral points in momentum space around 
which the dispersion is linear, leading to the so-called 
Dirac cones [l[. Much of the interest in graphene stems 
from an analogy of this dispersion relation to that of rela- 
tivistic, massless fermions, leading to intriguing new phe- 
nomena. For instance, a very unusual half-integer quan- 
tum Hall effect (QHE) and a non-zero Berry's phase have 
been discovered in graphene 0, 0, 0| , as well as an ab- 
normally weak localization 0,l6[. More recently, Raman 
spectroscopy 0, 0, 9| . flol| and angle resolved photoelec- 
tron spectroscopy 11( have been applied to graphene, 
yielding information on electron-phonon coupling and on 
the energy dispersion of the Dirac cones. 

Infrared (IR) spectroscopy is a powerful tool for inves- 
tigating the low-lying energy excitations of a material. 
When combined with a magnetic field, B, it allows for 
the study of its Landau level (LL) spectrum. In tra- 
ditional 2D materials with parabolic dispersions this is 
tantamount to measuring the carrier effective mass, to* , 
since transitions between the equally spaced LLs at en- 
ergy E n — (n + l/2)heB/m* reflect the same m* as in 
classical cyclotron resonance at lu c = eB/m*. Here e 
is the electron charge, Ti is Planck's constant, and the 
non-negative integer n is the LL index. 

In a magnetic field, the linear dispersion relation of 
graphene leads to an unequally spaced LL spectrum [l], 

mm, 



E n = sgn(n) x y / 2ehc 2 B \n\ = sgn(n) x y/2 \n\hc/lo (1) 



where c is the band velocity, Iq = ^Jh/eB is the magnetic 
length, and n > or n < represents electrons or holes, 
respectively. Most unusually, for n — there exists a LL 
at Eq = with a distinctive electron-hole degeneracy. 
This peculiar behavior of carriers in graphene is further 
enriched by spin splittings [H , possible lifting of the Dirac 
cone valley degeneracy and general many-particle effects, 
and provides a unique opportunity to probe these excep- 
tional electronic properties of graphene via LL formation. 

In this letter, we report IR transmission results on sin- 
gle layer graphene. In fields up to B = 18 T, two identifi- 
able LL transitions are clearly resolved, and their energy 
position scales as y/~B with a slope corresponding to a 
c « 1.1 X 10 6 m/s in Eq. (1). A deviation from an ideal 
ratio of 1 : (\/2 + 1) between the two transition ener- 
gies, predicted by simple matrix elements [LJL indicates 
a considerable many-particle contribution |l5[ to these 
LL transitions. 

Sadowski et al reported LL spectroscopy of ultrathin 
graphite layers created by thermal deposition on SiC [l6| . 
Their data - at considerably lower magnetic field - also 
show \f~B behavior and they deduce a c similar to ours. 
However, the nature of the thin graphite sheet remains 
unclear. 

Our studies were performed on graphene, mechani- 
cally cleaved from bulk Kish graphite and deposited onto 
lightly doped Si/Si02 substrates [l7| which are transpar- 
ent to IR, yet sufficiently conductive to serve as gates. 
The primary challenge of the experiment is the mismatch 
between the size of the IR focus in typical light-pipe 
systems (~1 mm) and the lateral dimension of typical 
mechanically extracted graphene samples (~10 fim). To 
overcome this limitation, we selected only large graphene 
specimens with areas as big as a few thousand /jm 2 and 
employed a parabolic cone to focus the IR light to a few 
hundred jim spot. Standard e-beam lithography, metal 
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FIG. 1: (color online). Half-integer quantum Hall effect in 
graphene with plateaus at R xy = h/ue 2 for filling factors 
v = ±2, ±6, ±10, at 4.2 K and 18 T. LL filling factors are 
indicated. Note that our devices are density tunable with 
n s /V g = 7.2 x 10 10 cm _2 V _1 , where n s is the charge carrier 
density in graphene. Upper inset: optical image of an 1100 
/im 2 graphene device. Lower inset: schematic of the experi- 
mental setup. 



evaporation and lift-off techniques were used to define 
the Cr/Au (3/35 nm) contact wires, along with a ~100 
/im diameter metal aperture around the specimen, which 
reduces stray IR light around the graphene sample, see 
inset to Fig. 1. The graphene device was mounted on the 
parabolic cone so the IR light passes through the Si sub- 
strate onto the graphene. The substrate was thinned and 
wedged to 4° to suppress Fabry-Perot interferences and 
could be aligned in situ with respect to the IR focus. 
A composite Si bolometer directly beneath the sample 
served as the IR detector. All electrical and IR experi- 
ments were performed at 4.2 K in magnetic fields up to 
B = 18 T using conventional quasi-dc lock-in technique 
and a Bruker IFS 66v/S FTIR spectrometer. 

The devices used in our work exhibit mobilities of 
2 — 4x 10 3 cm 2 V _1 s _1 measured at densities of ~ 2x 10 12 
cm~ 2 . The gate voltage corresponding to charge neu- 
trality, Voirac, is ^29 V due to a built-in potential from 
residual charges in the environment. Figure 1 shows a 
typical Hall resistance trace (Rxy) vs. gate voltage for 
fixed B = 18 T taken in situ. The characteristic half- 
integer QHE with plateaus at R xy — hjve 2 for filling 
factors v = ±2, ±6, ±10 0, 0| is clearly observed, con- 
firming the single layer nature of our graphene specimen. 

We record IR transmission spectra at fixed B field at 
two different carrier concentrations, corresponding to two 
different integer LL fillings. In this way the experimen- 
tal conditions are identical for all system components, 
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FIG. 2: (color online). Normalized IR absorption spectra 
of holes in graphene at three different magnetic fields, taken 
by dividing spectra taken at filling factors v = — 2 and v = 
— 10. Two LL resonances are denoted by Ti and Ta. Residual 
spectral artifacts are associated with 60 Hz harmonics and 
with carriers in the Si substrate. Dashed purple lines are 
Lorentzian fits to the data. The inset shows a schematic LL 
ladder with allowed transitions indicated by arrows. 



such as bolometer sensitivity, silicon transparency and 
B-induced shifts of optical path. The only change is in 
the LL occupation of the graphene sample. Furthermore, 
since the LL spacing - and hence the IR resonance - is 
density (LL filling factor) dependent, the two spectra, al- 
though at identical B-ficld, show transmission minima at 
two different IR frequencies and thus their ratio results in 
a maximum and a minimum placed on a background of 1. 
Figure 2 shows select transmission spectra of the sample 
shown in Fig. 1 at different magnetic fields following this 
normalization method, for two traces at filling factors 
v = —2 (numerator) and v = — 10 (denominator). Given 
the fourfold LL degeneracy of graphene and the electron- 
hole symmetry of the n — LL, filling factors v = —2 
and i/ — —10 correspond to a Fermi level (Ep) position 
between the n = — 1 and n = LLs and the n = — 3 
and n = —2 LLs, respectively; see inset to Fig. 2. At 
B = 18 T, switching the gate voltages between these two 
positions tunes the hole density from 9.4 x 10 11 cm -2 to 
4.7 x 10 12 cm" 2 . 

Two transmission minima, Ti and T2, are readily ob- 
servable in the traces of Fig. 2. Their minima associate 
them with the Ep{y = —2) position. Features from the 
Ep(y = —10) position occur as maxima residing at lower 
energies, and are only visible as shoulders in the 18T and 
12. IT traces. Therefore, in Fig. 2 the Ep(y = —10) spec- 
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tra simply serve as a normalization for the Ep(y = —2) 
spectra. 

Remaining artifacts in the data are largely due to 60 
Hz harmonics (narrow spikes). In addition, a feature 
associated with the Si substrate emerges, since the intro- 
duction of carriers in graphene leads to the same de nsity 
but opposite sign of carriers at the Si-SiC^ interface [l8( . 
The transmission of this carrier system is density depen- 
dent and therefore does not completely average out in 
spectral division. However, these artifacts in Fig. 2 can 
be circumvented and we find that the minima can be fit 
well by a Lorentzian, shown by the dashed lines. From 
these fits we determine the resonance position and the 
half-width of each resonance. 

The transmission minima in Fig. 2 are well developed 
and their resonance energies clearly decrease with de- 
creasing magnetic field. Evidently, T2 is considerably 
weaker than Ti and the ratio of their intensities seems 
to be roughly constant. The widths of the resonances, 
5E, are similar for Ti and T2 and correspond to a scat- 
tering time, t = K/5E = 20 fs in reasonable agreement 
with r = 40 fs, derived from the mobility. 

Figure 3(a) summarizes the energies of the Ti and T2 
transitions determined from Lorentzian fits to spectra 
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FIG. 3: (color online), (a) Resonance energies vs. vB, 
from holes (ratio of v = —2 and v — — 10 data, Fig. 2) 
and electrons (ratio of v = +2 and v = +10, spectra not 
shown) . Similar behavior has been observed in a second sam- 
ple (squares) . Solid line is a best vB-fit to the Ti transition, 
yielding £ = (1.12 ± 0.02) x 10 6 m/s. The dashed line rep- 
resents a scaling of the solid line by a factor (y2 + 1). (b) 
Half-width at half maximum from Lorentzian fits of the Ti 
transitions, as a function of \[B. 



from holes (Fig. 2) and electrons (not shown), plotted 
as a function of y/B. A clear \J~B relationship emerges. 
This dependence is very reliable for the Ti transition 
with a large number of data points to low i?-ficld. Ac- 
cording to Eq. (1) a s/B dependence is expected for any 
transition between LLs in graphene and our data con- 
firm this. Inspecting the slope and comparing it with 
Eq. (1) assuming c « 1 x 10 6 m/s @, lE E3, S iden- 
tifies the transition as the n = —1 — * n = (holes) 
and n = — » n = 1 (electrons) intraband LL transitions. 
Quantitatively, this assignment leads to a refined value of 
c = (1.12 ± 0.02) x 10 6 m/s, indicated as a line through 
the data in Fig. 3(a). Given the Fermi level position 
Ep(v = —2) shown in the inset to Fig. 2 these are the 
lowest LL transitions possible in graphene. Fig. 3(b) 
shows the -B-dependence of the half-width at half maxi- 
mum of transition Ti. It seems to be increasing, but its 
independence cannot be further assessed at this stage. 

Had this been a traditional 2D electron or hole system 
we would have observed only one resonance line, the cy- 
clotron resonance at Tiuj c = TteB/m*, from which arises 
a unique mass, to*. In contrast, we observe a second, 
higher energy transition, T2. To identify the LLs in- 
volved requires a recapitulation of the relevant dipole 
transition selection rules. In a traditional 2D system 
with energy spectrum E n = (n + l/2)hw c they dictate 
An = 712 —rii — ±1 for either electron or hole Landau 
ladders. When expanded to semiconductors in which in- 
tcrband transitions from a hole LL, rih, to an electron 
LL, n e , are also feasible, the selection rule changes to 
An = n e — n% = due to the separate p and s char- 
acter of valence and conduction band. In graphene, in 
which "valence" and "conduction band" have the same 
symmetry, a single selection rule An = [77.3 1 — |ni| = ±1 
is obtained [3] for intra as well as for interband tran- 
sitions. Therefore, the next allowed LL transitions for 
Ep(y — —2) with energy above the n = —1 — > n = 
transition are interband transitions n = — 2 — ► n = 1 
and n = — 1 — ► n = 2, see inset Fig. 2. While pre- 
serving the general y/B dependence their energy scales 
as 1 : (V2 + 1) compared to the n = — 1 — > n = intra- 
band transition, following Eq. (1). The dashed line in 
Fig. 3(a) represents a scaling of the Ti transition by a 
factor of y/2+ 1 and falls quite close to the T2 data. This 
provides very good evidence that T2 represents the low- 
est interband transition expected from the inset of Fig. 
2 and its electron- hole symmetric equivalent. 

The observed small deviation of the 1 : (v2+l) ratio of 
the Ti and T2 energies, however, lies well outside of the 
experimental errors (^symbol size) . A v^B-fit to T2 with 
a zero energy intercept results inc= (1.18 ± 0.02) x 10 6 
m/s, according to Eq. (1). This value differs appreciably 
from the value c = (1.12 ± 0.02) x 10 6 m/s deduced from 
Ti, as well as from the value of c = (1.03 ± 0.01) x 10 6 
m/s found by Sadowski et al [l6| . who derive the same 
band velocity c for both transitions. 
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The discrepancy between the band velocities deduced 
from different LL transitions sheds doubt on the appli- 
cability of a simple LL energy subtraction scheme based 
on Eq. (1) for the interpretation of IR data in graphene. 
In 2D systems with parabolic dispersions, Kohn's theo- 
rem [2lj explains that e-e interactions have no impact on 
the LL transition energies observed in IR experiments 
(cyclotron resonance). Instead, the resonance energy 
coincides with the non-interacting value, provided that 
the system is translationally invariant, which, in spite of 
residual disorder, can largely be assumed to hold. Yet, 
Kohn's theorem fails in the case of graphene, whose lin- 
ear dispersion may be viewed as a case of extreme non- 
parabolicity, and thus many-particle effects may be ex- 
pected to contribute to the LL transition energies. 

Indeed, the first calculations of many-particle correc- 
tions to the bare LL transitions are appearin g in the lit- 
erature. A recent paper by Iyengar et al [15[ arrives 
at quite large many-particle contributions, on the order 
of a few e 2 /elo = 56 meV \J B(T) /e, which amounts to 
an energy of ^60 meV at B = 18 T with an assumed 
dielectric constant of e = 4. Importantly, the total tran- 
sition energy between LLs with indices n and m scales as 
AE niTn xl = \/2xhc(^\m\±^y\n\)+C n ^ m xe 2 /e, so that 
contributions from c and C n ^ n j t cannot be distinguished 
by their independence. 

On the other hand, the magnitude of the many-particle 
effects, C n _ rn , depends on the specific n, m LL pair, a fact 
that can be employed to distinguish them from the band 
velocity contribution c. Iyengar et al [l5, 22 1 calculated 
the prefactors C_i.o = 1.18 and C_2 i = C_i 2 = 3.17 
for the transitions Ti and T2 of Fig. 3(a), respectively. 
Since C_2,i/C-i,o > (y/2 + 1), any e-e contribution to 
the LL transition will increase the ratio of transition en- 
ergies beyond its non- interacting value of V2 + 1, as is 
observed in our data. In fact, a band velocity c~8xl0 5 
m/s and e « 5 can fit simultaneously the Ti and T2 data 
of Fig. 3(a). However, at this stage a wide range of c 
and e produce similar good fits within the error bars of 
the experiment making it premature to assign any par- 
ticular values. More precise IR measurements between 
several other, higher LLs are required to assess quan- 
titatively the impact of many-particle physics on these 
transitions. However, it appears likely at this stage that 
a determination of the band velocity from a single parti- 
cle LL picture leads to a considerable overestimation of 
its value. A quantitative comparison of our data with 
many-body calculations [l5, 22 1 may be premature, since 
theory does not take into account disorders [23, 24|nor 
any mesoscopic corrugations of the graphene sheet [5j. 

Finally, apart from many-particle effects, one may 
wonder as to the impact of a possible gap opening around 
the Dirac point due to interaction effects, a scenario that 

r 1 

has been addressed by Gusynin et al [14| . In this model, 
any positive value for the gap energy reduces the ratio of 
T2 to Ti below v / 2 + 1, m contrast to our data. 



In summary, we have observed intra- and inter-LL 
transitions in IR spectroscopy on graphene. The tran- 
sition energies scale as \f~B~ and a simple, non-interacting 
LL transition interpretation yields a band velocity of 
c « 1.1 x 10 6 m/s. The observed deviation from a pre- 
cise 1 : (v / 2 + 1) scaling between the transition energies 
indicates a contribution from many-particle interactions 
to the transitions. Theory predicts rather strong (^30%) 
such corrections [H, HH . Qualitatively, the observed de- 
viations show the correct sign, but the magnitude of the 
corrections cannot yet be deduced reliably from experi- 
ment. 

We would like to thank Kun Yang, H.A. Fertig, L. 
Brey, Jianhui Wang, and D.N. Basov for discussions and 
B. Oezyilmaz for assistance in fabrication. This work 
is supported by the DOE (DE-AIO2-04ER46133 and 
DE-FG02-05ER46215), NSF (DMR-03-52738 and CHE- 
0117752), NYSTAR, and the Keck Foundation. IR mea- 
surement of this work was performed at the National 
High Magnetic Field Laboratory, which is supported by 
NSF Cooperative Agreement No. DMR-0084173, by the 
State of Florida, and by the DOE. We thank J. Jaroszyn- 
ski, E.C. Palm, T.P. Murphy, S.W. Tozer, and B.L. 
Brandt for experimental assistance. Z. Jiang and E.A. 
Henriksen contributed equally to this work. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 
[12] 
[13] 
[14] 

[15] 
[16] 
[17] 

[18] 
[19] 

[20] 



[21] 
[22] 
[23] 

[24] 



Electronic address: jiang@magnet.fsu.edu 
J.W. McClure, Phys. Rev. 104, 666 (1956). 
K.S. Novoselov et al., Nature 438, 197 (2005). 
Y. Zhang et al., Nature 438, 201 (2005). 
Y. Zhang et al., Phys. Rev. Lett. 96, 136806 (2006). 
S.V. Morozov et al., Phys. Rev. Lett. 97, 016801 (2006). 
Xiaosong Wu et al., Phys. Rev. Lett. 98, 136801 (2007). 
A.C. Ferrari et al., Phys. Rev. Let t. 97, 187401 (2006). 
D. Graf et al |cond-mat/0607562| 
J. Yan et al., |cond-m at/0612634 



aL,| cond-mat/0611714] 



Simone Pisana et 

Aaron Bostwick et al., cond-mat/0609660 

F. D.M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

G. W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 
V.P. Gusynin, S.G. Sharapov, and J. P. Carbotte, 
|cond-mat/0607727 



A. Iyengar et aL,|cond-mat /0608364I 
M.L. Sadowski et al., Phys. Rev. Lett. 97, 266405 (2006). 
K.S. Novoselov et al., Proc. Natl. Acad. Sci. U.S.A. 102, 
10451 (2005). 

N. Sai et al., Phys. Rev. B 75, 045307 (2007). 

M.S. Dresselhaus and G. Dresselhaus, Adv. Phys. 51, 1 

(2002). 

N.B. Brandt, S.M. Chudinov, and Y.G. Ponomarev, 
Semimetals 1: Graphite and Its Compounds (North- 
Holland, Amsterdam, 1988). 
W. Kohn, Phys. Rev. 123, 1242 (1961). 
Jianhui Wang and H.A. Fertig, private communication. 
N.M.R. Peres, F. Guinea, and A.H.C. Neto, Phys. Rev. 
B 73, 125411 (2006). 

Balazs Dora and Peter Thalmeier, cond-mat/0701714 



